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We discuss results of 2D simulations of magnetorotational(MR) mechanism of core collapse supernova explosions. Due 
to the nonuniform collapse the collapsed core rotates differentially. In the presence of initial poloidal magnetic field 
its toroidal component appears and grows with time. Increased magnetic pressure leads to foramtion of compression 
wave which moves outwards. It transforms into the fast MHD shock wave (supernova shock wave). The shape of the 
MR supernova explosion qualitatively depends on the configuration of the initial magnetic field. For the dipole-like 
initial magnetic field the supernova explosion develops mainly along rotational axis forming mildly coUimated jet. 
Quadrupole-like initial magnetic field leads to the explosion developing mainly along equatorial plane. Magnetorota- 
tional instability was found in our simulations. The supernova explosion energy is growing with increase of the initial 
core mass and rotational energy of the core, and corresponds to the observational data. 
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1. Introduction 

The problem of explanation of the core-collapse su- 
pernova event is one of the interesting and long 
standing problems in astrophysics. Mechanisms, 
based on the bounce shock energy or neutrino in- 
teraction with the matter of pre-supernova star do 
not lead to the supernova explosion. 

The MR mechanism for core collapse super- 
nova explosion was suggested by Bisnovatyi-Kogan 
in 1970.""^ The main idea of the MR mechanism is 
to transform part of the rotational energy of pre- 
supernova into the radial kinetic energy (explosion 
energy). Due to the not uniform collapse the iron 
core rotates differentially. Differential rotation leads 
to the appearing and amplification of the toroidal 
component of the magnetic field. Growth of the mag- 
netic field means amplification of the magnetic pres- 
sure with time. A compression wave appears near the 
region of the extremum of the magnetic field. This 
compression wave moves outwards along steeply de- 
creasing density profile. In a short time it transforms 
to the fast MHD shock wave. When the shock reaches 
the surface of the presupernova it ejects part of the 
matter and energy of the presupernova star. This 
ejection can be interpreted as an explosion of the 
core collapse supernova star. First ID simulations of 
MR supernova explosion were represented in,^ see 
also.^ 

The first simulations of the MR processes in stars 



have been done by,^ after which MR processes in the 
stars in relation to the core collapse supernova explo- 
sion had been simulated by.^'^"^ Recently the inter- 
est to the MR processes (especially in application 
to the core collapse supernova) was recommenced 
due to increasing number of observational data about 
asymmetry of the explosion, and possible coUimated 
ejecta in connection with cosmic gamma ray bursts 

^9-14,16^ 

Our results of simulations of MR supernova ex- 
plosion mechanism show that this mechanism allows 
to produce 0.5 — 1.3 ■ lO^^ergs energy of explosion. 
These values of SN explosion energy correspond to 
estimations made from core collapse SN observa- 
tions. 

It was found, that the shape of the MR SN ex- 
plosion qualitatively depends on the configuration 
(symmetry type) of magnetic field. The initial field 
of quadrupole type of symmetry leads to the MR 
explosion which develops predominantly near equa- 
torial plane, while the initial dipole type field results 
in SN explosion as mildly coUimated jet developing 
along axis of rotation. 

MR instability was revealed in 2D simulations of 
MR SN mechanism. MR instability leads to the ex- 
ponential growth of both poloidal and toroidal com- 
ponents of the magnetic field and significantly reduce 
time of the MR supernova explosion in comparison 
with ID simulations.^ 

For the 2D MHD simulations we used specially 
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Fig. 2. Example of velocity field for the MR supernova ex- 
plosion with the initial quadrupole field. 



Fig. 4. Example of specific angular momentum distribution 
for the MR supernova explosion with the initial dipole field. 



developed numerical code based on the implicit com- 
pletely conservative Lagrangian scheme on triangu- 
lar grid of variable structure (see"'^^ and references 
therein). It was tested thoroughly with different 
tests. The example of the Lagrangean triangular gird 
is given at the Fig.l. 

Grid reconstruction procedure applied here for 
the reconstruction of the triangular lagrangian grid 
is used both for the correction of the "quality" of the 
grid and for the dynamical adaptation of the grid. 



2. Formulation of the problem 

For the simulations of the magnetorotational super- 
nova with initial dipole-likc magnetic field leading to 
the formation of protojet, we used the set of MHD 
equations with self-gravitation and infinite conduc- 
tivity as in the papers. ^^'^^ We rewrite the MHD 
equations in nondimensional form using the same 
scales as in the papers. ^^'^^ We used the same equa- 
tion of state and neutrino losses formulae as in the 
papers, and have started calculations from the 
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same model as in the papers. ^^'^^ 

The initial poloidal magnetic field is defined as 
in our previous paper^^ by the toroidal current 
using Bio-Savara law. The toroidal current which de- 
termines the initial magnetic field slioidd be defined 
in the upper and in the lower hemispheres. The field 
with the quadrupole-like symmetry is formed by the 
toroidal current antisymmetrical to the equatorial 
plane. The dipolc-like magnetic field is formed by 
current symmetrical to the equatorial plane. 

2.1. Basic equations 

Consider a set of magnetohydrodynamical equations 
with self-gravitation and infinite conductivity: 

dx 



dt 



dv ^/ H-H 

d /H 



+ 



dp 
dt 

V • (H (g) H) 



4-77 



pV • V = 0, 



H • Vv, A$ = AwGp, 



hprp -i — 
"^^^ dt ~ dt 



X : 



p^+PV-v + pF{p,T) = 0, 
P = P{p,T), e = e{p,T). 
V • V is the total time derivative, 



[Vr,V, 



Vz) is the velocity vector, p 



is the density, P is the pressure, H = {Hr, H^, Hz) 
is the magnetic field vector, $ is the gravitational 
potential, e is the internal energy, G is gravitational 
constant, H (g) H is the tensor of rank 2, and F{p, T) 
is the rate of neutrino losses. 

r, v'j and z are spatial Lagrangian coordinates, 
i.e. r = r(ro,</?o, and zo,i), V = viro.Vo.zo.t), and 
z = z{rQ,^fQ,ZQ,t), where ro, <^0) zq are the initial co- 
ordinates of material points of the matter. 

Taking into acc;onnt symmetry assumptions 
=0), the divergency of the tensor He 
presented in the following form: 



) H can be 



V-(H(g)H) 



1 djrH^H^) 

r dr 

1 O(rHrH^) 

r dr 
r dr 



+ 



dz 
dz 
dz 



Axial symmetry = 0) and symmetry to the 
equatorial plane are assumed. The problem is solved 
in the restricted domain. At t = the domain is 
restricted by the rotational axis r > 0, equatorial 
plane z > 0, and the outer boundary of the star 



where the density of the matter is zero, while poloidal 
components of the magnetic field H^, and Hz can be 
non-zero. 

At the rotational axis (r = 0) the following 
boundary conditions are defined: (V$)r = 0, Vr = 0. 
At the equatorial plane {z = 0) the boundary condi- 
tions are: (V$)z = 0, Vz = 0. At the outer boundary 
(boundary with vacuum) the following condition is 

defined, ^outer boundary — 0* 

3. MR supernova explosion 

MR explosion with the initial quadrupole-like mag- 
netic field was described in detail in the paper. Af- 
ter the core collapse the pre SN rotates differentially. 
The toroidal component of the magnetic field ap- 
pears and grows linearly with time at the initial stage 
of the MR explosion. When the toroidal magnetic 
field reaches some certain value its linear growth 
changes to the exponential growth of the toroidal 
and poloidal components due to the appearing of the 
magnetorotational instability (MRI). A toy model 
for the qualitative explanation of the MRI in MR 
supernova was suggested in.^^ MR supernova explo- 
sion with initial quadrupole-like magnetic field re- 
sults in explosion which develops mainly near equato- 
rial plane. At the Fig. 2 the velocity field and specific 
angular momentum Fig. 3 for the initial quadrupole- 
like field. 

Simulations of the MR supernova with the ini- 
tial magnetic field of dipole-like symmetry leads to 
the qualitatively different result in the shape of ex- 
plosion.^'' In this case the MR explosion develops 
mainly along the axis of rotation and forms mildly 
collimated proto jet (Figure 4). The protojet found 
in our simulations could be collimated when it de- 
velops in the extended envelope of the massive star 
- the progenitor of the core collapse supernova. 

MRI leads to the formation of chaotic magnetic 
field structure. I the case of finite conductivity the 
reconnection of the magnetic field could be impor- 
tant. We have estimated characteristic time of the 
reconnection of the magnetic field using results of 
our simulations.^^ We found that characteristic time 
for the reconnection of the magnetic field for super- 
nova parameters used in our simulations is « 5s. The 
MR supernova explosion time in our simulations is 
w 0.5— Is what is significantly less then characteristic 
time of the magnetic field reconnection development, 
and do not influence the MR supernova significantly. 
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Fig. 5. Dependence of the supernova explosion energy on the 
core mass for initial angular velocity uJa ~ 3.53s^^ (solid line) 
and uiQ ~ 2. 525"^^ (dashed line) (before collapse) 

Wc have done simulations for different initial 
masses of the iron core and for different initial angu- 
lar velocities of presupernova. The supernova explo- 
sion energy grows significantly with increase of the 
core mass. The dependence of the explosion energy 
on the core mass for the different initial values of the 
rotational energy (angular velocity) is presented at 
Fig. 5 for the initial quadrupole-like magnetic field. 
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4. Conclusions 

The rcst.ilts of 2D simulations of MR supernova ex- 
plosion mechanism show that it allows to get the 
explosion energy which corresponds to observational 
values of the explosion energy for core collapse su- 
pernovae. The magnetorotational instability appear- 
ing in the simulations of MR supernova significantly 
reduces the explosion time, and increase the chaotic 
magnetic field in the young neutron stars. "'^^ Wc have 
found that MR supernova explosion energy is grow- 
ing with increase of the initial mass of the iron core 
and initial rotational energy. 
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